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. We discuss an approach for studying the properties of mesoscopic systems, where discrete and 

' continuum parts of the spectrum are equally important. The approach can be applied (i) to stable 

\ heavy nuclei and complex atoms near the continuum threshold, (ii) to nuclei far from the region of 

nuclear stability, both of the regions being of great current interest, and (iii) to mesoscopic devices 
with interacting electrons. The goal is to develop a new consistent version of the continuum shell 
model that simultaneously takes into account strong interaction between fermions and coupling to 
the continuum. Main attention is paid to the formation of compound resonances, their statistical 
properties, and correlations of the cross sections. We study the Ericson fluctuations of overlapping 
resonances and show that the continuum shell model nicely describes universal properties of the 
conductance fluctuations. 
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(~| , A. Introduction 

O . The focus of interest in low-energy nuclear physics has recently concentrated on nuclei far from the region of 
stability. The success of the nuclear shell model with strong interaction between nucleons urges one to look for 
ways of incorporating this rich experience into a more general context that would properly include the continuum. 
' The continuum shell model is expected to combine the physics of intrinsic structure and reactions in contrast to 
the traditional spectroscopy and reaction theory which were historically developed as different branches of nuclear 
I . physics. In this approach, the properties of open channels, threshold behavior, symmetries, and spectroscopic factors 
■ O ' enter the overall theory, so that a fully consistent picture uniting structure and reaction mechanisms is created. With 
great successes of mesoscopic physics (complex atoms and molecules, micro- and nano-devices of condensed matter, 
atomic and molecular traps, and prototypes of quantum computers), the nuclear problem becomes a typical part of 
general science of open mesoscopic systems. 

In this paper we discuss a new approach [l| for studying collective phenomena in nuclei and quantum dots which 
emerge as a result of strong interactions between fermions coupled to the continuum. In this problem the standard shell 
model techniques, adjusted for the bound states, cannot work and continuum effects must be properly accounted for. 
At a certain strength of the continuum coupling, the system of strongly interacting fermions undergoes a collective 
^\ restructuring. This is an area where the conventional division of nuclear physics into "structure" and "reactions" 
y—i becomes inappropriate and two views of the process, from the inside (structure and properties of bound states) and 
' , from the outside (cross sections of reactions), have to be unified. 
' Conventional theory of statistical reactions does not answer the question of interplay between reactions and internal 
00 ' structure determined by the character of interactions between the constituents. Here more detailed considerations are 
required based on the generalization of the shell model of nuclear reactions [2] . Such an extension requires statistical 
^ [ assumptions concerning intrinsic dynamics and the coupling to the continuum [sj, To account for specificity of 
. the system, one has to go beyond standard random matrix approaches 0, 0, [71 utilizing the Gaussian Orthogonal 
' Ensemble (GOE). The consistent description based on the continuum shell model 3| , as well as more phenomenological 
approaches [^, [^ , indicates the presence of a sharp restructuring of the system when the widths of resonances become 
comparable to their energy spacings. This phenomenon carries features of a quantum phase transition with the 
strength of continuum coupling playing the role of a control parameter. 

As was numerically observed in the shell model framework [1(J|, the distribution of resonance widths rapidly changes 
in the transitional region in such a way that a number of very broad resonances equal to a number of open decay 
channels absorb the lion's share of the total width of all overlapped resonances, while the remaining states become 
very narrow. The corresponding theory was suggested in [Til . Il2l . [isj , where the mechanism of this restructuring 
was understood to be associated with the nature of the effective non-Hermitian Hamiltonian [2] that describes the 
intrinsic dynamics after eliminating the channel variables. The segre gati on of short-lived broad resonances from long- 



lived trapped states was shown to be similar to the superradiance in quantum optics induced by the coupling 
of atomic radiators through the common radiation field, an analog of coherent coupling of many overlapped intrinsic 
states through continuum decay channels. Later a general character of the phenomenon was demonstrated for systems 
with the GOE intrinsic dynamics and many open channels [isl . [l6j . Modern versions of the shell model in continuum 
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[ITI . Il8j are based on the effective Hamiltonian and naturally reveal the superradiance phenomenon as an important 
element. The transition to this regime should be taken into account in all cases when a physical system is strongly 
coupled to the continuum, see for example [l^ and references therein. 

In the present work, using the same framework as in we study the interplay between the intrinsic dynamics 
and statistical properties of cross sections comparing the results with those of conventional approaches, namely 
Hauser-Feshbach average cross sections and Ericson fluctuations and correlations. We also analyze how the universal 
conductance fluctuations appear in our model, how they depend on the intrinsic chaos and what is the role of 
correlations between cross sections of various processes. 



B. The model 



We consider a generic system of n interacting fermions distributed over m single-particle levels of energies ( "mean- 
field orbitals"). The fermions are coupled to continuum through the finite number M of open decay channels. In 
the many-body representation, the Hamiltonian matrix of size N = m!/[n!(m — n)!] entirely determines the intrinsic 
dynamics of the closed system, and our main interest is to reveal how the properties of the reactions depend on 
the de gree of intrinsic chaos that is due to the inter-particle interaction. According to the well developed formalism 
[3, IT^. Il9l|. the behavior of an open system can be described in terms of the effective non-Hermitian Hamiltonian 
consisting of two terms, 



M 

Here if is a Hermitian matrix that can be also influenced by the presence of the continuum [13, [H, H^l, and W 
stands for the coupling between intrinsic many-body states labeled as and open decay channels a, 6, c... The 

factorized structure of W is dictated by the requirement of the unitarity of the scattering matrix. In what follows we 
assume that the transition amplitudes A\ between intrinsic states \i) and channels c are real quantities, therefore, we 
restrict ourself by the time-invariant systems. As a result, both H and W are real symmetric matrices. 

The Hermitian part, H = iJo + of the full Hamiltonian ([T]) is modeled by the so-called two-body random 
interaction (TBRI) (see, for example, [2l[ and references therein), 

m ^ 

H = Ha + V; i/o = X! ^ ^ i^'^ysis.^s^.s^al.al^as^.as^. (2) 

s=l 

Here the mean- field part Hq describes non-interacting particles (or quasi-particles), and the interaction between 
particles is given by V. In the many-particle basis \k) the matrix H is constructed by the Slater determinant, 
= a\_^ . . . a^^ |0), where a\ and as are the creation and annihilation operators. The interaction between the 
particles is considered to be of two-body nature, therefore, each many-body matrix element Vik = {l\V\k) is a sum of 
a number of two-body matrix elements VsiS2S3S4 involving at most four single-particle states \s). In our simulations 
we take n = 6, m = 12 that provides a sufficiently large dimension N = 924. 

The single-particle energies, Cs, are assumed to have a Poissonian distribution of spacings, with the mean level 
density l/d^. The interaction V is characterized by the variance of the normally distributed two-body random matrix 
elements, (ysi,s2\s3,si) — ''^o- Without the interaction, wq = 0, the many-body states have also the Poissonian spacing 
distribution P{s). In the opposite limit of an extremely strong interaction, X = vq/cIq oo (or, equivalently, for 
do = 0), the function P(s) is close to the Wigner- Dyson (WD) distribution that is typical for completely chaotic 
systems Q . The critical value of the interaction for the onset of strong chaos was obtained in Ref. f2l| , it is estimated 
as 

_ Vcr _ 2(m - n) 
da Ns 

where iV, = n{m — n) -I- n{n — 1) (m — n) (m — n — l)/4 is the number of directly coupled many-body states in any row 
of the matrix Uij. One can expect that for a very strong interaction, A ^ Acr, some of properties of the system may 
be close to those obtained with iJ taken from the Gaussian Orthogonal Ensemble (GOE). For this reason, we have 
also considered, for comparison, the GOE matrices in place of H. 

In our study the amplitudes A\ are assumed to be random independent Gaussian variables with the zero mean and 
variance 

(A^A/)=<5,,<5-'^. (4) 
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This is compatible with the GOE or TBRI models where generic many-body states coupled to continuum have a very 
complicated structure, while the decay probes specific simple components of these states related to few open channels. 
Even for a weak intrinsic interaction, one needs to have in mind that in reality the states \i) have certain values of 
exact constants of motion, such as angular momentum and isospin in the nuclear case. Therefore, at sufhciently large 
dimension N, these states acquire geometric chaoticity [2^ due to the almost random coupling of individual spins. 
Thus, the assumption of a random nature of the decay amplitudes is reasonable. 

The parameters 7"^ characterize the total coupling of all states to the channel c. The normalization used in Eq. 
(HI is convenient if the energy interval ND covered by decaying states is finite. Here D is the distance between the 
many-body states in the middle of the spectrum, D — l/p{0), where p{E) is the many-body level density, and E = 
corresponds to the center of the spectrum. Below we neglect a possible energy dependence of the amplitudes that is 
important near thresholds and is taken into account in realistic shell model calculations [l^, . The ratio ^'^ /ND 
characterizes the degree of overlap of the resonances in the channel c. The control parameter determining the strength 
of the coupling to the continuum can be written as follows, 

(5) 
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The transition from separated to strongly overlapped resonances corresponds to k"^ w 1. In order to keep the coupling 
to continuum fixed, in our numerical calculations we renormalize the absolute magnitude of the widths, 7"^, by varying 
the intrinsic interaction and, therefore, the level density p. 

The effective Hamiltonian allows one to construct the scattering matrix, 

gab ^ ^ab _ -^ab^ (g) 

where 



r^\E)=^An^-^) A] (7) 

are the scattering amplitudes determining the cross sections cr°'^{E) of reactions, 

a^^E) = \T"'\E)\'^. (8) 

In our notations the cross sections are dimensionless since we omit the common factor i: /h?' . In what follows we study 
both the elastic, b = a, and inelastic, b a, cross sections. Note that we ignore the smooth potential phases that are 
irrelevant for our purposes. 



C. Transition to superradiance 

The complex eigenvalues £r — (*/2)rr of Ti. coincide with the poles of the S'-matrix and, for small 7"^, determine 
energies and widths of separated resonances. In the simulations we consider a real energy interval at the center of the 
spectrum of Ti. with the constant many-body level density p(0) = D^^. The transmission coefficient in the channel c, 

T^ = l-|(5->P = -i^, (9) 

is maximal (equal to 1) at the critical point, k"^ = 1, and, for simplicity, we assume M equiprobable channels, k'^ = k. 
For each value of k we have used a large number of realizations of the Hamiltonian matrices, with further averaging 
over energy. 

As 7'^ grows, the resonances start to overlap, leading to a segregation of the resonance widths [HI [H, HI] 
occurring at k « 1, see Fig. 1. The widths of M resonances are increasing at the expense of the remaining N — M 
resonances. For a weak coupling, k <C 1, the widths are given by diagonal matrix elements, Fi = (i = J^^ii^i)^ , 
and the mean width is (F) = ^M/N. In the limit of strong coupling, k ^ 1, the widths of M broad resonances converge 
to the non-zero eigenvalues of the matrix W that has a rank M due to its factorized structure. As for the remaining 
"trapped" (iV — M) states, their widths decrease as I/7. From our data one can conclude that the transition to the 
superradiance in dependence on the control parameter k turns out to be quite sharp, and may be treated as a kind 
of the phase transition. 
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FIG. 1: Distribution of widths, F, and positions, uj, of resonances (poles of the scattering matrix) as a function of k. The results 
are shown for n = 6 particles, m = 12 orbitals, 50 open channels, and for a relatively strong interaction between particles, 
Vol do = 1/10, above the strong chaos threshold. The smooth curves are the boundaries due to the analytical results obtained 
in Ref. [2^ for the GOE matrices in place of H 
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FIG. 2: Normalized average width versus k for n = 6, m = 12 and M = 10. The average was performed over the resonances 
around the center of the energy spectrum. Open circles refer to the GOE case, pluses to A = 0, squares to A — > cx), and crosses 
to A = 1/30. Solid curve shows the MS-expression (|10|l . Full circles stand for the normalized correlation length at A = 1/30, 
see in the text. The dashed curve shows the Weisskopf relation (fT6)) . 



In Fig. 2 we show how the average width (normalized to the mean level spacing D) depends on k for different values 
of the interaction A. The figure demonstrates a sharp change in the width distribution at the transition point k = 1. 
The results are compared with the expression 
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(10) 



that is obtained from the Moldauer-Simonius (MS) relation [1, H^l for M equivalent channels. We have extrapolated 
the MS-relation to the region with k > 1 for which only N — M narrow resonances were taken into account. For one 
channel in the GOE case this result was obtained in Refs. [HI, [13) [iBl ; here we see that the MS-expression is also valid 
for a large number of channels. Note that the results are independent of the degree of intrinsic chaoticity; this fact 
was explained in Ref.fH. According to [l^, the divergence of (F) at k = 1 is due to the (non-integrable) power-law 
behavior for large F, see below; in the numerical simulation, see Fig. 2, there is no divergence because of the finite 
number of resonances, although the trend is clearly seen. 

As for the width distribution P{T), in Ref. [l| it was shown that for a weak coupling, k ^ 1, the conventional xif- 
distribution is valid for any strength of the interaction. However, when k grows, a clear dependence on the interaction 
strength emerges. As first noted in Ref. [25], as k increases, P{T) becomes broader than the Xm distribution. For 
large k, both for the GOE and for A — > cx), P(r) is again given by the Xai distribution, contrar y to the cases of the 
finite interaction strength. Our results are in a good agreement with the analytical predictions jl6l . [26j obtained for 
the GOE case (neglecting the emergence of broad resonances in the regime of strong overlap at k > 1). In particular, 
for large number of channels, the tails of P(r) decrease as F~^, thus leading to the divergence of (F) at k = 1. 
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D. Ericson fluctuations 

According to the conventional theory of Ericson fluctuations [27l . [28l . [29j , the scattering amphtude can be written 
as the sum of the average and fluctuating parts, T''''{E) = {T''^{E)) + T^''{E), with {T^''{E)) = 0. Correspondingly, 
the average cross section, a = |Tp, can be also divided into two contributions, (a) = ((Tdir) + (ffl)- Here the direct 
reaction cross section, (crdir), is determined by the average scattering amplitude only, while (aa) is the fluctuational 
cross section (also known as the compound nucleus cross section). 

The main assumption of the standard theory is that in the regime of overlapping resonances, (F) > D, the fluctuating 
amplitudes can be written as T^'^{E) ~ £, + irj where both ^ and rj are Gaussian random variables with zero mean. 
This is supported by the fact that for (F) 3> D both ^ and rj are the sums of a large number of random variables. This 
assumption can be associated with the other one, namely, with statistical independence of poles (resonance energies) 
and residues (resonance amplitudes). As a result, the statistical methods can be developed that lead to various 
conclusions. In the frame of our model we have tested the most important of the Ericson predictions concerning the 
resonance widths and cross sections (see also Refs.[l|), paying special attention to the dependence on the intrinsic 
interaction strength, A. 

a. Fluctuations of resonance widths In the theory of Ericson fluctuations it is usually assumed that the fluctua- 
tions of resonance widths are small for a large number of channels, w{r)/{T)'^ <ti 1, where w(r) stands for the variance 
of widths. In justiflcation of this statement [27j it is argued that in the overlapping regime the width of a resonance 
can be presented as a sum of M partial widths. Assuming that they obey the Porter-Thomas distribution, the total 
width is expected to have a xti distribution, so that w{T)/ (F)^ = 2/M is small for M ^ 1. 

However, recently we have found fl*] that for large values of k the distribution of the widths strongly differs from the 
xii distribution. The data show that, as k increases, the normalized variance, ti;(F)/(F)^, also increases, remaining 
very large even for M = 20 1. Thus, contrary to the traditional belief, the relative variance of widths does not 
become small for large number of channels. The data confirm that for non-overlapping resonances (small k <^ 1) 
there is no dependence on A and w{r)/ (F)^ decreases as 2/AI for all the ensembles, as expected. However, as k grows, 
the dependence on A emerges: the weaker the intrinsic chaos (and, consequently, the more ordered is the intrinsic 
spectrum) the larger are the width fluctuations. For large k ^ 1 and strong interaction between particles, A > Acr, 
the width distribution is again given by the Xai distribution, contrary to the cases of the finite interaction strength. 

b. Cross sections Speaking about cross sections and their fluctuations, it is convenient to use the so-called elastic 
enhancement factor, that for M equivalent channels takes the form, 

where h ^ a. Here {a^) stands for the fluctuational part of cross sections (elastic for a — b and inelastic for a ^ b). 
In the case of equal channels, we obtain 

,ab. ^ i~\{sn\' ^ T 

^ ^ ' F + M -I F + M -V ^ ' 

therefore, ((Tg°) = F{a^^), where T is the transmission coefficient, see Eq. Since the transmission coefficient T 
does not depend on A, the only dependence on A in Eq. ([12]) is contained in the elastic enhancement factor F. The 
same seems to be correct even when the channels are non-equivalent, according to the results of Ref. [s^]- It should 
be pointed out that F also depends on k. Specifically, with an increase of k from zero, the value of F decreases, being 
confined by the interval between 3 and 2, see discussion in Ref. fT]. 

As one can see from Ea.([T2|). with an increase of the number of channels the dependence on the interaction 
strength A disappears for the fluctuational inelastic cross section. In contrast, the fluctuational elastic cross section 
manifests a clear dependence on A, according to the data of Ref. This fact allows one to directly relate the value 
of the enhancement factor F to the strength A of interaction between the particles. Specifically, the more regular is 
the intrinsic motion, the higher is the average cross section. For a large number of channels, the A-dependence of the 
elastic cross section is in agreement with the estimate, {<t^°') FT/M . 

c. Fluctuations of cross sections The fluctuations of both elastic and inelastic cross sections strongly depend on 
the coupling to the continuum, see examples in Fig. 3. According to the standard Ericson theory, in the region of 
strongly overlapping resonances, k w 1, the variance of fluctuations of both elastic and inelastic cross sections can be 
expressed via the average cross sections. Specifically, for equal channels one gets. 



Var(a) = {an) (2(adir) + {(Jr)) = {(Jti)^- 



(13) 
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FIG. 3; Elastic cross-section versus energy E for different strength of coupling to the continuum: for weak, k = 0.07, strong, 
K — 0.7, and very strong, k — 7.0, for M = 10. Other parameters are the same as in Fig.l. 



The latter equality is written for the absence of direct processes, that is our case. Thus, one can expect the rela- 
tion, Var((T)/((Tfl)^ — 1, for a large number of channels. Our data confirm this relation for inelastic cross sections, 
independently of the interaction strength A, and for M > 10. 

As for the elastic cross sections, the expression ([13]) emerges for a quite large number of channels, M > 25. This 
fact of slow convergence to the limiting value was also found analytically in Refs. [sU, [s^l for the GOE case, any 
number of channels and any coupling strength with the continuum. However, simple expressions were derived only 
for MT ^ 1. Our data reveal a clear dependence on A that is stronger for the non- normalized variance Var(cr) than 
for the normalized variance Var(cr)/ (crfl)^. This effect was explained in Ref. [l|. 

Another well known prediction of Ericson theory is that for strongly overlapping resonances both inelastic and 
elastic cross sections have the exponential distribution, 

P(x)=e-^ (14) 

Indeed, our data for inelastic cross sections show good agreement with this prediction, independently of the interaction 
strength A and for large number of channels, M > 10. For the elastic cross sections, the exponential distribution arises 
for larger number of channels, M > 25, with a weak dependence on A. In general, one can say that the exponential 
form of the fluctuations is a good approximation in the region k w 1. 

d. Correlation functions The correlation function of cross sections is defined as 

Cie) = {aiEME + e))-{aiE)r. (15) 

By assuming the Gaussian form of distribution for fluctuating amplitudes 7fl, one can obtain that the correlation 
function has a Lorentzian form, with the correlation length equal to the average width, la- — (F). In our model we 
have found that l^ ~ Is for a large number of channels in the region k w 1, and for any interaction strength A. Here, Is 
is the correlation length obtained from the correlation function of the scattering matrix. On the contrary, for smaller 
M, our data show that Is < la, and this difference grows for the weaker interaction between the particles. 

For large M, the correlation functions have, indeed, the Lorentzian form for all A. Meanwhile, for a small number 
of channels the correlation function is not Lorentzian, in agreement with the results of [s^. It is important to stress 
that for any M, the correlation length is different from the average width, apart from the region of small k. Instead, 
for a large number of channels, the correlation length is determined by the Weisskopf relation, see jsj and references 
therein, 

/ _ MT _ M 4k 

Z) ~ ^ ~ 2^(H-k)2- ^^^^ 

The Weisskopf relation is confirmed by our data in Fig. 2, it has been also derived in Ref. fs'] for small values 
of the ratio m = M/N , in the overlapping regime for the TBRI model with the infinite interaction, as well as in Ref. 
Q for the GOE ensemble. The fact that the correlation length is not equal to the average width was recognized long 
ago, see Refs. [11] and Q. 



E. Conductance fluctuations 



Our model can be also used to study conductance and its fluctuations for quantum dots with interacting electrons. 
In such an application one can treat M/2 channels as the left channels corresponding to incoming electron waves, and 



other M/2 channels, as the right channels corresponding to outgoing waves. Then, one can define the conductance G 
in a standard way, 

M/2 M 

G = E E ^'^^ (17) 

a=l 6=Af/2+l 
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FIG. 4: Average conductance versus «; for tlie GOE and very strong interaction, A oo (left panel), and for a weak interaction, 
A = 0.02 (right panel) for different number of channels, M — 4, 10, 16. Dashed curves correspond to the analytical expression 
P8|) with F = 2.0, (left panel), and dotted curves with F = 2.5, (right panel). The data are obtained for n = 6 particles, 
m = 12 orbitals, and M — 4, 10, 16 channels. 

Therefore, the properties of the conductance are entirely determined by the inelastic cross-sections, h ^ a. In our 
case of equivalent channels the average conductance reads, 

(G)^—{'T )-—F + M-i^— ^^^^ 

where T = 4k/ (1 + k)^ and the last expression is valid for M ^ 1. As one can see, the influence of the inter-particle 
interaction is due to the enhancement factor F only. This factor ranges from F — 3 for the absence of the interaction, 
A = 0, to F = 2 for a very strong interaction that can be modeled by the GOE ensemble [sBl- Since currently there 
is no theory relating the enhancement factor to the degree of chaos in the internal dynamics, we used this factor as 
the fitting parameter. In Fig. 4 we report our numerical data obtained for the average conductance in dependence 
on the coupling strength k for different values of M. As expected, for a very strong interaction (absence of the mean 
field part, do = 0, or, the same, A oo) the results are very close to those obtained for the GOE Hamiltonian H, 
and there is a good agreement with the expression ^Tl} with the theoretical value F — 2. On the other hand, for a 
relatively weak interaction, A = 0.02, (below the chaos threshold in a closed model) the data can be described by the 
same expression with F = 2.5 that is in the middle of the region between the Poisson and Wigner-Dyson statistics 
for the energy spacings. As one can see, for this intermediate case the correspondence between the relation and 
data is also good. Note, that for a very large number of channels the influence of the internal dynamics disappears. 

Our main interest was in the variance, Var(G), of the conductance when we change the energy E around the center 
of the energy spectrum. In Fig. 5 we report the data for the variance for the GOE case (left part) and for a relatively 
weak interaction A — 0.02 in dependence on the coupling parameter k and number M of channels. First observation 
concerns a clear evidence that for a weak coupling, k 0.1, the role of the interaction should not be neglected. This 
fact is in contrast with the average conductance for which the influence of the internal chaos seems to be much less 
in comparison with the conductance fluctuations. 

Second important observation is that at the transitional point, k « 1 (known in the solid state applications as 
"perfect coupling"), practically there is no influence of the inter-particle interaction. This fact is quite instructive, 
it manifests an important role of the superradiance when the (strong) interaction between the particles is mainly 
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FIG. 5: Variance of the conductance versus n for the GOE (left panel), and for A = 0.02, (right panel) . Circles, squares and 
crosses stand for M — 16, 10, 4 channels, respectively. The data are given for n — 6 and m = 12. 



due to the continuum. Also, the dependence of fluctuations on the number of channels at the transitional region is 
relatively weak (and disappears with an increase of M). In the theory of conductance fluctuations it is known that for 
a perfect coupling the variance takes the universal value in the limit of a very large number of channels. This value 
for time-invariant systems is either 2/15 or 1/8 depending on the type of transport (ballistic or diffusive, both for a 
bulk scattering, see, e.g. [131 )■ It is quite instructive that our data confirm the value 1/8 for the GOE at k = 1 with 
a very good accuracy. Specifically, the data also reveal small corrections due to the finite values of M, that reduce 
the variance to that of slightly below 1/8, in agreement with the analytical results obtained in Ref. [33|)- Also, the 
data indicate the value 2/15 for the TBRI model with a moderate interaction between the particles. 

We have performed a detailed analysis of the emergence of the universal value 1/8. We should point out that one 
has to take into account the correlations between different cross sections in the expression for the variance of the 
conductance. Otherwise, one gets the incorrect result. 



Var(G) = — Var(a) = — 



T 



F + M -1 



4 ' 



(19) 



where the last expression is obtained for AI ^ 1. Thus, for the perfect coupling, we would have Var(G) ~ 1/4, which 
is twice the correct value. 

In order to reveal the role of correlations between different cross sections, instead of Eq. (fT9)) we write the relation 
that takes into account all the correlations. 



/ T 

Var(G) = — = 

^ ' 4 \F + M -1 



+ N*C^ + {N,-N*)Ca, 



(20) 



where Nc — L{L — 1), N* = L{M — 2), L — M^/4. Here the terms Ca and Gs stand for the correlation functions, 

G = (o-«''o-«'''') - ((T«'')(cr'''^') (21) 

with one common index (either a = a' or 6 = b') for the E— correlations, and with no common indexes (both a ^ a! 
and h ^ h') for the A— correlations. The bar in Eq. ((2T|) represents the average over different possible combinations 
of a, 6, a', 5', and (...) indicate the average over energy. Our analysis shows that for M ^ 1 one obtains, Gs « —M~^ 
and Ga « 2M~^. Therefore, in the limit of large number of channels, one can write, in the correspondence with the 
structure of Eq.([20l). 



Var(G) = i-i + i.i. 
^ ' 4 4 8 8 



(22) 



This result clearly demonstrates the crucial role of correlations for the variance of the conductance. These correlations 
are neglected in the traditional theory of Ericson fluctuations. 

Of special interest is the k— dependence of the correlations for a non-perfect coupling, in the regions k < 1 and 
K > 1. As one can see, the A— and S— correlations have very different behavior, both for weak and strong intrinsic 
chaos. The remarkable fact is that, for the perfect coupling with k — Q, the A— and E— correlations have maximal 
values of opposite signs. These results may be important for the analysis of experimental nuclear or solid state data. 
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FIG. 6: Dependence of average correlations between cross sections on the coupling n for the GOE case (left panel) and for a 
weak interaction with A = 0.02 (right panel). Circles stand for the Ce and for the Ca correlations (refer to text for details). 
Both panels are on the same scale for n = G and m = 12. 

F. Conclusions 

In conclusion, we have studied the statistics of resonance widths and cross sections for a fermion system coupled 
to open decay channels. We used the effective non-Hermitian Hamiltonian that reflects exact quantum-mechanical 
features of many-body dynamics in an open interacting system. Main attention was paid to various signatures of the 
crossover from isolated to overlapping resonances in dependence on the strength of inter-particle interaction modeled 
by the two-body random matrices in the intrinsic part of the non-Hermitian Hamiltonian. 

In the framework of our model we have tested main predictions of the Ericson fluctuation theory. We have shown 
that the assumption that the fluctuations of the resonance widths become negligible for a large number of channels, is 
invalid in the overlapping regime. We found that the fluctuations of resonance widths increase along with the coupling 
to the continuum, resulting in the divergence of the relative fluctuation of the widths (the ratio of the variance to the 
square of the average width) at the transition point k = 1 for any number of channels. 

We have also shown that, in agreement with previous studies, the correlation length differs from the average 
width for any number of channels. On the other hand, the Weisskopf relation that connects the correlation 
length of the cross section to the transmission coefficient, works, for a large number of channels, at any value of the 
intrinsic interaction strength A. In many situations the data show that increase of A suppresses the fluctuations in 
the continuum. This can be understood qualitatively as a manifestation of many-body chaos that makes all internal 
states uniformly mixed [38]. 

Our results can be applied to any many- fermion system coupled to the continuum of open decay channels. The nat- 
ural applications first of all should cover neutron resonances in nuclei, where rich statistical material was accumulated 
but the transitional region from isolated to overlapped resonances was not studied in detail. Other open mesoscopic 
systems, for example, quantum dots and quantum wires, should be analyzed as well in the crossover region. 

It was our special interest to study the mesoscopic fluctuations of the conductance in the dependence on the coupling 
strength and degree of intrinsic chaos. In particular, we have analyzed how the correlations between different cross 
sections influence the variance of the conductance fluctuations. Our data manifest that these correlations determine 
the universal value 1/8 (see, for example, (stI. [39j) for the conductance variance in the perfect coupling regime, k — 1, 
and for very strong interaction between particles. Thus, we demonstrate that the Ericson theory of cross-section 
fluctuation that neglects these correlations, is insufflcient for the description of the mesoscopic fluctuations. 
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